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We consider the average number B m (t) of bonds traversed exactly m times by a i step simple 
random walk. We determine B m (t) explicitly in the scaling limit t — > oo with m/\/t fixed in 
dimension d — 1 and m/logt fixed in dimension d — 2. The scaling function is an erfc in d — 1 
and an exponential in d — 2. 
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I. INTRODUCTION 

The simple random walk is a venerable problem, finding its applications in many areas of statistical physics. Despite 
its long history, novel aspects continue to surface. Motivated by recent experiments on gas transport properties in 
polycarbonate films, in which polymers are found in a glassy state, we are led to the the following question. Suppose 
the polymers are modeled by the usual (non-interacting) random walks on a lattice, how many "monomers" can 
we expect on each bond? Clearly, the answer will depend on both the length and the density of the polymers. Not 
surprisingly, it is a simple step once we know the distribution of bond-traverses of a single polymer (or random walker) . 
Remarkably, this distribution is not, to the best of our knowledge, in the literature. This note will be devoted to our 
findings of such a distribution. 

Let a simple random walk start at the origin of a d-dimensional hypercubic lattice. We wish to know the average 
number of bonds B m (t) that it has traversed exactly m times. This quantity satisfies 

oo oo 

Y J B rn {t)=B{t), ^mB m {t)=t (1) 

m—1 m—1 

where Bit) is the average total number of distinct bonds traversed by (also known as the support) the walk. 

If for t — ► oo asymptotically B(t) ~ (3(t), then one may reasonably expect that in that limit B m (t) can be expressed 
as a scaling function of m(3(t)/t. The sum rules (|lj) then imply that B m (t) takes the form 

where B is a scaling function. 

In this note we will explicit find B (and j3) in spatial dimensions d = 1 and d = 2, thereby justifying the scaling form 
Eq. (|^). It appears that B is an error function "erfc" in d = 1 and an exponential in d = 2. In d = 2 we determine, 
moreover, the leading order correction to the scaling behavior (|^). 



II. GENERATING FUNCTION METHOD FOR REPEATED BOND TRAVERSALS 



In order to determine B m (t) wc need the following notation. Let the vector 5 denote any of the d basis vectors of 
the lattice. Let (a;, S) denote the bond between the sites x and x + 5 (the fact that this bond is identical to (se + 5, —5) 
is of no importance). Let B m (x,S;t) be the probability that (x, S) is traversed exactly m times, irrespective of the 
direction. Then 
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B m (t) = J2 B m(x,S;t) (3) 

(x.S) 

where each bond occurs exactly once in the summation. The quantity B m (x, 8] 4) may be calculated by an adaptation 
of the analogous method for multiple visits to the same site Q. Let F(x,S;t) be the probability that the first traversal 
of (x, 8) occurs at the 4 th step, with 4 = 1,2,... It is convenient to set F(x, 8; 0) = 0. 

For any X(t) defined for 4 = 0,1,2,... we may introduce the generating function X(z) = z* X (t) . Let 

furthermore R(t) be the probability that, given a traversal has taken place, the next one occurs exactly r steps later, 
with r = 1,2,... Then R(t) = F(0,8;t) and hence R(z) = F(0,6;z). Reasoning in a similar way as for the site 
problem one finds that 

B m (z) = ^— z [E P ^ *5 z )] t 1 ~ k ^ ^-'(z) (m = 1, 2, . . .) (4) 

Here the function F, which implies R, still have to be found. 

Let G(x; 4) be the probability that after 4 steps the walk is at lattice site x, for 4 = 0, 1, 2, . . . Let G(x, 5; t) denote 
the probability that at its ith step it traverses (in either direction) the bond (x, 5), for t = 1,2,... We set additionally 
G(x, 6; 0) = 0. Since (x, S) can be traversed starting either from x or from x + S, we have 

G(x,6;t) = ^[G(x;t-l)+G(x + S;t-l)] (4=1,2,...) (5) 

By a slight extension of the standard procedure for calculating first passage probabilities on sites |^,[l],|| we have here 

t 

G(x,S;t)=F(x,S;t) +^2F(x,S\T)G(0,S]t-T) (4 = 1,2,...) (6) 

In terms of generating functions, Eqs. (||) and (||) become, respectively, G(x, 5; z) — (z/2d)[G(x; z) + G(x + 5; z)] and 
F(x, 5; z) — G(x, 5; z)/[l + G(0, 5; z)]. Elimination of G(x, 6; z) from this pair of equations yields 

p ± G(x;z) + G(x + S;z) 

M 1 + ^[G(0;z) + G(S;z)] 

This achieves the reduction of the desired function F(x,S;z) to the known function G(x;z). Note that expression 
(0) has the required invariance under the replacement (x, 8) — > (x + 8, —8). We now exploit the well-known relations 
G(8; z) = [G(0; z) — l]/z and J2 X G{x\ z) = 1/(1 — z). Upon combining these with Eqs. ([|) and (^) one gets, writing 
henceforth G(z) = G(0; z), a fully explicit expression for the generating function B m (z) in terms of G{z). 



(2d) 2 r zG(z) 
(1 - zf [2d + zG(z)} 2 \-2d + zG(z) 



B m {z) = n hr 2 7Z7^k ~, n : 'I ; . (8) 



where G (z) is defined by 



The m-th coefficient is extracted as 



G = z- 1 [{l + z)G(z)-l] 



1 f d 

B m {t) = —j-^- l B m {z), (9) 

where the integral runs counterclockwise around the origin. To make B m (t) more explicit, we have to consider each 
spatial dimension separately. 
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III. ONE DIMENSION 



In dimension d = 1 we have to evaluate Eq. (||) with ||,[T| G(z) = (1 — z 2 ) 1 / 2 . It turns out to be advantageous to 
consider the differences B m+ i(t) — B m (t). After slight rewriting this yields 



B 



rn+1 



(t) - B m {t) 



1 

7ri 



dz 



yt+m+2 



i-vT 



(10) 



We will now fold the integration path around the branch cut that runs from z = 1 to z = oo along the positive real 
axis. In the limit t — > oo we expect a meaningful result only if also m — * oo at fixed ratio m/y/i. We anticipate that 
in this limit the integral on z draws its main contribution from a region at a distance of order i -1 / 2 from the branch 
point at z = 1. We therefore introduce the scaling variable 



/Lti = m/V2t (11) 

and the scaled variable of integration y = (z— l)t. The integrand of Eq. (|l^) may now be expanded in powers of t^ 1 / 2 
at fixed /j and y. This leads to 



B n 



,(t)-B m {t) 




y/V 



(12) 



where the two terms on the right hand side come from above and below the branch cut, respectively. The integral 
is easily found to be equal to — 2 3 / 2 (7rt) -1 / 2 e _Ml . With the boundary condition B^t) = we therefore find upon 
integrating 

B m {t) =2erfc(/ii) + O^ 1 / 2 ) (13) 

which is the final result, valid for t — > oo at \i\ fixed. When summing Eq. (fi~3) on m (or alternatively when evaluating 
the integral obtained by summing Eq. (^) on m), one obtains B(t) ~ \/&bpK. Together with Eq. ( |l3[ ) this confirms the 
validity of the hypothesized scaling form (||). 

In order to compare the large time scaling behavior of B m (t) to its finite time forms we have performed Monte-Carlo 
simulations. For t = 10 2 , 10 3 , 10 4 the bond distribution was averaged over several independent runs. Fig. 1 shows 
that for increasing values of t the simulation data rapidly converge to the scaling function, and they are practically 
indistinguishable on this figure for t > 10 3 . 



pq 




m / 1 



FIG. 1. Simulation results for B m (t) in one dimension at different finite values of time and its scaling form in the t — + oo 
limit (solid line). 



IV. TWO DIMENSIONS 



For a two-dimensional square lattice the expansion G(z) = — log y^-j + C((l — z)log(l — z)) enables one to 
evaluate Eq. (||) for t asymptotically large. The only integral whose asymptotic behavior we need is j|]3) 
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1 

2?ri 



yt + l 



1 - Z f log™ -J 



ir n t 
log™ 8t 



n (1 - C) 
log 8i 



O(log- 2 80 



(14) 



where C = 0.577215... is Euler's constant. We anticipate a meaningful result for i — > oo if also m — ► oo in such a way 
that the ratio to/ logi remains fixed. We expect that in this scaling limit the z integral will draw its main contribution 
from a region in the complex plane near the branch point z = 1 where the ratio m/G(z) is finite. Keeping this in 
mind we expand the integral (0) for B m {t) as 



B m [t) 



1 
27ri 



dz 



s t+i (i 



z) 2 G 2 (z 



„-2m/G( Z ) 



1 



/TO \ 



G(z) 



+ 0(G- 2 (z)) 



(15) 



We may now insert in ( ^5[ ) the explicit expression ~ ^\ogj^, knowing that the 0((l — z)log(l — z)) terms 
in the expansion of G{z) will contribute only terms with higher powers of t^ 1 to the final asymptotic series. The 
exponential in Eq. (|l5|) may then be expanded, the resulting series integrated term by term with the aid of Eq. Jl^), 
and summed again. Upon introducing the scaling variable 



/'2 



2nm 
log8i 



one finds the final result 



B m {t) = 



4n 2 t 
log 2 8t 



- + 1 

1 + (M2 - 2) 2 + 



C 



log8t 



(log" 2 8t) 



(16) 



(17) 



valid for t — > oo at fixed. As expected, the correction terms decay only as powers of the inverse logarithm of 
the number of steps. The next few higher order terms in the asymptotic series ( p"7| ) may be calculated without great 
effort. Since the fcth order correction term is multiplied by a fcth degree polynomial in /12, this is actually an expansion 
in powers of /X2/log8t. One may verify that Eq. (|l7|) satisfies the sum rules ([!]) up to and including the first order 
correction term. The scaling function B introduced in Eq. ^ here appears to be a simple exponential. One finds that 
the average total number of bonds traversed increases as (3(t) = 27rf/log8i + 0{t log- 2 8t). 




FIG. 2. Simulation data for B m (t) in two dimensions at different values of time (symbols) compared to the theoretical scaling 
form (solid lines). The exact result for t — 10 2 (dashed line) and the scaling law plus leading correction for t — 10 8 (dotted 
line) are also displayed. 

Fig. 2 compares simulation data for t = 10 2 , 10 3 , . . . , 10 8 to the theoretical scaling law (the first term of Eq. (fL7|)). 
The t = 10 2 data are also compared to the exact curve, which we obtained by Taylor expanding Eq. (||) through 
the one hundredth term with the aid of a symbol manipulation program. For t — 10 s the scaling law plus leading 
correction (the first two terms of Eq. (|l7|) ) is also displayed. Upon collapsing the data of Fig. 2 one finds Fig. 3, 
which shows that in d = 2 the convergence to the asymptotic scaling law is very slow. 



4 





10 15 
27tm / log(8t) 



FIG. 3. Data collapse of two dimensional simulation results for B m (t). The t — > oo scaling limit is shown as a solid line. 



V. CONCLUDING REMARKS 



We conclude this note with a series of comments and remarks. 

An alternative route to B m (t) is to consider a complementary process, by focusing on a particular bond in, say, a 
finite periodic L d lattice. Now, we can ask how often this bond is traversed by a t-step random walk which starts 
at xq and ends at x. By summing over (x,Xq) and invoking translational invariance, we obviously access B rn {t). In 
this approach, we can easily generalize the problem to a study of the frequency of traversing in only one direction. 
To solve both problems, consider a modified random walker, for which the rate of traversing our chosen bond is 
p/2d instead of l/2d. Solving the master equation for P(x,t\xo,0;p) (the usual probability) by standard generating 
function techniques and defining B (p, z) = J2 X x St^o zt P{ x i A x 0i 0;? 5 ); we obtain 



B(p,z) 



L" 



2z(l-p) 



1 - z 



(1 - zf 2d + (1 - p) zG (z) 



(18) 



The first term represents, since ~}2, x P(x,t\xQ,§;p = 1) = 1, the L d points of origination of the walks we considered. 
Due to this difference in normalization, the bivariate generating function for B m (t), i.e. Em=i E™oP mzt -^ m Wi * s 
precisely d B (p, z) — B (0, z) . With the extra factor d, the quantity dB (1, z) represents nothing but the total number 
of bonds in the lattice. Interestingly, this B (p, z) contains an extra term (B (0, z)) which carries the information on 
the bonds never traversed. So, d L d / (1 — z) — B (0, z) is just zt B(t). For completeness, we report the result 

for unidirectional traverses: 



B(p,z) = 



L" 



z(l-p) 



(l-zy[2d+{l-p)zG(z)} 



Finally, note that finite size effects are fully incorporated in this approach, although they are implicitly "buried" in 

r , -i-l 

G (z) = L~ d J2{k } 1 — ( z /d) Y^i—i cos (where the sum is over the set {ki} of allowed L d wavevectors) . Needless 
to say, generalizations to strip- or slab-like samples (L\ < < ...) are straightforward. Thus, it is possible to study 
the crossover of our distributions, at least in principle, when the polymer length exceeds the shortest dimension 
significantly (t 3> L\). Physically, ultra-thin membranes made from extra- long polymers can be manufactured. Until 
they become reality, however, it may not be worthwhile to extend our analysis to this class of crossover behavior. 

A second remark concerns the statistics of multiple visits to sites, as opposed to traversals of bonds. Visits have 
been of considerable interest in the literature (see Hughes jj] and references cited therein). The average number V m (t) 
of sites visited exactly m times by a t step random walk was studied by Montroll and Weiss |5| and by Barber and 
Ninham B. The t — > oo limit at fixed m was considered by Hughes Q| (see also 0) in d = 1 and by Montroll and 
Weiss [[| (see also |8)) in d = 2; however, the scaling limit expression of V m (t) has not to our knowledge appeared in 
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the literature. Since the analysis of V m (t) runs exactly parallel to that of B m (t), we content ourselves to state the 
results here. 

In dimension d = 1 the average number of sites visited m times is, to leading order in the scaling limit, equal to 
the average number of bonds visited m times: V m (t) — 2erfc(/ii) for t — > oo at [i\ fixed. In dimension d = 2 one has 
in terms of the scaling variable /j.' 2 — 7rm/log8i 



7T t 



log^ 8t 



l + C 



v m(t) = ; — e " M2 1 + (M2 " 2) 2 = p , + 0(log~ 2 8i) (19) 



log 8t 



valid for t — > oo at ^ 2 fixed. Comparison of Eqs. ( |l9|) and ( |l7| ) shows that to leading order V m (t) and B m (t) are 
identical up to a coefficient and a scale factor; however, the coefficients of the first correction terms are different. The 
leading order relation between V m (t) and B m (t) in d = 1, 2 is heuristically clear as follows. Given a large number m 
of visits to an arbitrary site x, there will have been typically m/ (2d) traversals, starting from x, of a specific bond 
(x, x + 5). Now the same bond will have been traversed, typically, the same number of times in the opposite direction. 
So for each site visited m times, there are d bonds traversed m/d times. 

Although we have not shown so explicitly, the scaling function in d = 2 is expected to be universal, i.e. lattice 
structure independent. In dimensions d > 2 the random walk is transient and it is easy to show that as a consequence 
in the large t limit B m (t) ~ b m t and V m (t) ~ v m t , where b m and v m are nonuniversal. For random walks that are not 
simple (i.e. have a step size distribution not limited to nearest neighbor steps), bond traversals are not unambiguously 
defined. However, visits to sites still are, and for such walks V m (t) is readily calculated by the present method. One 
case of interest is the scaling limit of V m (t) for lattice Levy flights (also called Riemann walks) |§,[l(J, for which the 
distribution of step sizes decreases as a power law. 
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